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Abstract In the present article, we investigate the prop- 
erties of motion in Hamiltonian systems of two and 
three degrees of freedom, using the distribution of the 
values of two new dynamical parameters. The distribu- 
tion functions of the new parameters, define the S(g) 
and the S(w) dynamical spectra. The first spectrum 
definition, that is the S(g) spectrum, will be applied in 
a Hamiltonian system of two degrees of freedom (2D), 
while the S(w) dynamical spectrum will be deployed in 
a Hamiltonian system of three degrees of freedom (3D). 
Both Hamiltonian systems, describe a very interesting 
dynamical system which displays a large variety of res- 
onant orbits, different chaotic components and also sev- 
eral sticky regions. We test and prove the efficiency and 
the reliability of these new dynamical spectra, in detect- 
ing tiny ordered domains embedded in the chaotic sea, 
corresponding to complicated resonant orbits of higher 
multiplicity. The results of our extensive numerical cal- 
culations, suggest that both dynamical spectra are fast 
and reliable discriminants between different types of or- 
bits in Hamiltonian systems, while requiring very short 
computation time in order to provide solid and conclu- 
sive evidence regarding the nature of an orbit. Further- 
more, we establish numerical criteria in order to quan- 
tify the results obtained from our new dynamical spec- 
tra. A comparison to other previously used dynamical 
indicators, reveals the leading role of the new spectra. 
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1 Introduction 

Over the last half century, a large number of studies 
have been dedicated to efforts on finding new ways to 
characterize the nature of orbits in Hamiltonian sys- 
tems [1,2,7,18,22,26,30]. Among them the dynamical 
spectra (i.e. distribution of the values of a given param- 
eter along the orbit) have been proved a very reliable 
tool for the exploration of the properties of motion in 
dynamical systems. In previous works [37] and [11-12], 
presented a new method that allows the distinction be- 
tween domains of ordered and chaotic motion in Hamil- 
tonian systems of two and three degrees of freedom, or 
in 2D and 4D symplectic maps. This method is based 
on the calculation of stretching numbers, helicity an- 
gles and twist angles (see the above referenced papers 
for theoretical details) and their respective dynamical 
spectra. 

Consider the Poincare section map of a Hamiltonian 
dynamical system. Let be the next image of an in- 
finitesimal deviation ^ between two nearby orbits. The 
"stretching number" ai [37], or the "local Lyapunov 
indicator" [16] is defined as 



oii = In 



i+l 



(1) 



i.e. it is a 1-step Lyapunov Characteristic Number (LCN). 
The spectrum of stretching numbers is their distribu- 
tion function, namely 



S(a) 



AN(a) 
NAa 



(2) 



where AN(a) is the number of stretching numbers in 
the interval (a, a + Aa), after N iterations. Other types 
of dynamical spectra have also been used in previous 
studies [3,5,6,9,10,19,20,28]. 
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In the present paper we will demonstrate and prove 
the high efficiency of the new dynamical spectra in sepa- 
rating ordered from chaotic domains in a realistic galac- 
tic model. The separation can be obtained even with an 
integration time for individual orbits, as short as 5000 
time units. Thus, our new dynamical indicators are 
faster than other spectral methods, e.g. the frequency 
analysis method [21], and much more faster than the 
calculation of the Lyapunov Characteristic Exponent 
(LCE) [23]. In particular, our new dynamical spectra 
are very fast and practical tools for locating multiple 
islands of invariant curves embedded in the chaotic sea, 
which correspond to resonant orbits of higher multi- 
plicity. In addition, they can be used to provide a quick 
estimation of the fraction of the phase space volume 
occupied by ordered or chaotic orbits, at many dif- 
ferent energy levels. This information is needed, when 
one wishes to construct self-consistent models, in which 
many different orbits must be combined. 

We shall test our new dynamical spectra in the po- 
tential 



The Hamiltonian corresponds to the potential (3) is 



V to t{x,y,z) = V har {x,y,z) + V nuc (x,y,z), 



where 



V har (x, y,z) = ^- (x 2 + ay 2 + bz 2 ) 



and 



Vnuc(x> V: Z) 



M n 



^Jx 2 + y 2 + z 2 + c 2 n 



(3) 



(4) 



(5) 



Here a and b are the flattening parameters, M n and c„ 
is the mass and the scale length of the nucleus respec- 
tively, while the parameter oj is used for the consistency 
of the galactic units. We use a system of galactic units, 
where the unit of mass is 2.325 x 10 7 M Q , the unit of 
length is 1 kpc and the unit of time is 0.997748 x 10 8 
yr. The velocity unit is 10 km/s, while G is equal to 
unity. In the above units we use the values: u> = 10 km 
s^kpc" 1 , a = 4,6 = 1.25, M n = 50, while c n = 0.25. 
The 3D gravitational potential V tot , describes satisfac- 
torily the central parts of an elliptical galaxy and con- 
sists of two parts: the first part, that is V} lav is the 
potential of a 3D anisotropic harmonic oscillator, de- 
rived from Taylor expansion of the logarithmic poten- 
tial, while the second part, V nuc , is a Plummer poten- 
tial of a spherical nucleus. This particular dynamical 
system has been studied in detail in [38] (hereafter Pa- 
per I). The main reason which justifies the choice of 
the potential (3) is that it displays a large variety of 
resonant orbits, different chaotic components and also 
several sticky regions. Therefore, it provide us an ex- 
cellent opportunity in order to test and also reveal the 
great efficiency of our new dynamical parameters. 



H = \ (pI+pI+pI) + V to t(x,y,z) = h, 



(6) 



where p x ,p y and p z are the momenta per unit mass 
conjugate to x, y and z respectively, while h is the nu- 
merical value of the Hamiltonian. All the outcomes of 
the present research, are based on the numerical inte- 
gration of the equations of motion 



dV t 



x = — - 



tot 



dx 



u 2 + 



3/2 



dVto 
dy 



dV t 



tot 



dz 



bu 2 + 



(x 2 + y 2 + z 2 + c 2 n ) 

Mn 

(x 2 + y 2 + z 2 + c 2 n f 2 
M n 

(x 2 +y 2 + z 2 + c 2 n f 12 



*>(7) 



where the dot indicates derivative with respect to the 
time. We integrated the equations of motion with a dou- 
ble precision Runge-Kutta algorithm of 7 - 8 th order 
(RK78), in Fortran 77 [31]. The accuracy of the out- 
comes, was checked by constancy of the energy integral 
(6), which was conserved up to the fifteenth significant 
figure. 

Our aim is to test and explore the properties and 
also the advantages of our new dynamical spectra. The 
S(g) spectrum will be applied in the dynamical system 
of two degrees of freedom (2D), while for the dynamical 
system of three degrees of freedom (3D), described by 
the Hamiltonian (6), will shall use the S(w) spectrum. 
The layout of this article is as follows: In Section 2, we 
present an analysis of the phase plane of the 2D Hamil- 
tonian system. Furthermore, we apply the S(g) dynam- 
ical spectrum in order to study the various resonant or- 
bits and the different chaotic components. In the same 
Section, the evolution of the 2D sticky orbits is inves- 
tigated using the S(g) dynamical spectrum. In Section 
3, we investigate the properties of the 3D Hamiltonian 
system, using the new S(w) dynamical spectrum. Once 
more, we test and prove the efficiency of our new indi- 
cator by studying the 3D resonant orbits, the different 
chaotic components and also the sticky regions. Finally, 
we close with Section 4 in which we present a discussion 
and the conclusions of the present research. 

2 Application of the S(g) dynamical spectrum 
in the 2D Hamiltonian system 

In this Section, we shall study the nature of orbits in the 
2D dynamical model, using the S(g) spectrum. The cor- 
responding Hamiltonian for the 2D system (z = p z = 0) 
is 

H2 = Upl+p 2 y )+Vt ot (x,y) = h 2 , (8) 
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Fig. 1 The (x,p x ), y = 0, p y > phase plane, when = 1.5. 
Details are given in the text. 

where h 2 is the numerical value of 2D Hamiltonian. 

Figure 1 shows the (x,p x ), y — 0, p y > Poincare 
phase plane for the Hamiltonian (8) when h 2 = 1.5. 
The values of the involved parameters are: u> = 10, a = 
4, M n = 50 and c„ = 0.25. We observe, that the phase 
plane shown in Fig. 1 is very interesting and compli- 
cated. There is a plethora of sets of islands of invariant 
curves which correspond to resonant orbits of lower and 
higher multiplicity. Moreover, we discern the presence 
of several different chaotic components and also consid- 
erable sticky regions, embedded in the unified chaotic 
domain. These interesting features of the phase plane 
will help us to test, verify and eventually prove the ef- 
ficiency and the reliability of the new dynamical spec- 
trum. 

It is therefore evident from Fig. 1, that this dynam- 
ical system presents a large variety of resonant orbits, 
which correspond to sets of multiple islands of invari- 
ant curves in the (x,p x ) phase plane. Figure 2a-h shows 
eight representative periodic orbits for the 2D system. 
The values of all the other parameters are as in Fig. 1. 
Figure 2a shows a periodic orbit characteristic of the 
1:1 resonance with initial conditions x = 0.2943 and 
VxQ = 0. For all 2D orbits ?/o = 0, while the value of 
p y0 is always found from the energy integral (8). In 
Figure 2b one may observe a 1:2 periodic orbit with 
initial conditions xq = 0.548 and p x o = 10.8. Figure 2c 
depicts a 2:3 resonant periodic orbit with initial condi- 
tions x = and p x0 = 13.07. Moreover, in Figure 2d 
we can see a periodic orbit of the 3:4 resonance. The 
initial conditions are xo = 0.097 and p x o = 7.7. Figure 
2e presents a 3:5 resonant periodic orbit with initial 
conditions x — 0.889 and p x0 — 0. Similarly, in Figure 
2f a 4:5 resonant periodic orbit is presented. The initial 
conditions are xq = and p X Q = 4.675. In Figure 2g we 
observe a complicated resonant periodic orbit with ini- 



tial conditions xq = 0.667 and p X Q = 0, which belongs 
to the family of the 5:7 resonance. Finally, a periodic 
orbit of the 5:8 resonance is presented in Figure 2h. The 
initial conditions are x = 0.803 and p xa = 2.03. The 
integration time for all orbits shown in Fig. 2a-h is 100 
time units. It is remarkable, that all the above resonant 
families of orbits exist simultaneously in this 2D Hamil- 
tonian system. The interaction of all these resonances, 
justifies the presence of the different sticky regions and 
also the different chaotic components displayed in the 
phase plane of Fig. 1. 

In Figure 3a we present with different colors the 
most significant sets of islands of invariant curves in the 
(x,p x ) phase plane. These islands of invariant curves 
correspond to different types of families of resonant or- 
bits. The large central island is marked in black and 
corresponds to the 1:1 resonance. The initial conditions 
are xq = 0.15 and p X Q = 0. The islands of invariant 
curves corresponding to the 1:2 resonance are red and 
the initial conditions are xo = 0.5 and p x o — 12.25. The 
islands of invariant curves which correspond to the 2:3 
resonance are colored in green. The initial conditions 
are xo = and p x o = 14.6. Moreover, the 3:4 resonant 
islands with initial conditions xq = 0.1 and p X Q = 8.5 
arc plotted in blue. The set of islands of invariant curves 
which corresponds to the 3:5 resonance is marked in yel- 
low. The initial conditions are x — 0.909 and p x0 = 0. 
The islands of invariant curves produced by the 4:5 res- 
onance are depicted with magenta color and the initial 
conditions are x — and p x o = 4.021. The set of 
islands of invariant curves corresponding to 5:7 reso- 
nance is orange. The initial conditions are x = 0.657 
and p x o = 0. Last but not least, the islands of invari- 
ant curves produced by the 5:8 resonance are plotted 
in brown and the initial conditions are x = 0.802 and 
p x0 = 1.44. The outermost black solid line shown in 
Fig. 3a is the Zero Velocity Curve (ZVC). The integra- 
tion time for all the sets of invariant curves shown in 
Fig. 3a is 2000 time units. The three different chaotic 
components together with the unified chaotic sea are 
shown in Figure 3b. The first chaotic component (CC- 
1) has initial conditions x = —0.82 and p x0 — 0, the 
second (CC-2) x — —0.86 and p x o = 0, while the ini- 
tial conditions for the third chaotic component (CC-3) 
are Xq = — 1 and p x o = 0. Here, we notice that the 
third chaotic component coincides with the ZVC. All 
the chaotic components and the chaotic sea presented 
in Fig. 3b were integrated for 2000 time units. The equa- 
tion of the limiting curve - ZVC (that it the curve con- 
taining all the invariant curves of the 2D system) is 

\pl + V tot (x) = h 2 . (9) 
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Fig. 2 (a-h): Different types of resonant periodic orbits in the 2D dynamical system. The initial conditions and the values of 
all the involved parameters are given in the text. 
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Fig. 3 (a-b): (a-Ze/t): The most significant sets of islands of invariant curves in the (x,p x ) phase plane and (b-right): the three 
different chaotic components together with the unified chaotic sea. (Color figurer online) 



In an earlier work [39] (hereafter Paper II) we intro- 
duced and used a new dynamical parameter the S(g) 
spectrum, in order to study the islandic motion of res- 
onant orbits and the evolution of sticky regions. The 
S(g) spectrum was applied in a axially symmetric galac- 
tic model in the meridian (r, z) plane. However, it is 
easy enough to modify its definition in order to make 
the S(g) spectrum applicable in Cartesian coordinates 
(x, y). Thus, we define the dynamical parameter g as 

'i "T Pxi ^iPxi 



fji 



Pyi 



(10) 



where (xi,p X i,p y i) are the successive values of the (x,p x , 
elements of the 2D orbit on the Poincare (x,p x ), y = 0, 
p y > phase plane. We shall define the dynamical spec- 
trum of the parameter g, through its distribution func- 
tion 

AN {g) 



Py) 



S{g) = 



NAg 



(11) 



where AN{g) are the number of the parameters g in 
the interval (g, g + Ag), after N iterations. By defini- 
tion, the g parameter is based on a combination of the 
coordinates and the velocities of the 2D orbit. It is very 
interesting to test the effectiveness of the S (g) spectrum 
in the present dynamical system and find out if it has 
the ability to identify all the different types of resonant 
orbits. Moreover, it would be extremely beneficial, if it 
could distinguish the different chaotic components and 
also if it is sensitive enough in order to provide reliable 
results regarding the evolution of the 2D sticky orbits. 

First we shall apply the S(g) dynamical spectrum, 
in order to identify the different types of resonant or- 
bits which correspond to the sets of islands of invariant 
curves shown in Fig. 3a. The left side of Figure 4a-h 
shows the S(g) spectrums for the resonant orbits, while 
the right side of the same figure presents the corre- 
sponding islands of invariant curves in the (x,p x ) phase 
plane. Figure 4a shows the S(g) spectrum for the 1:1 



resonant orbit with initial conditions as in Fig. 3a. As 
expected, one can observe a well defined "t/-shaped" 
structure with two abrupt peaks in the two extreme op- 
posite values of g, indicating regular motion. The cor- 
responding island is shown in Figure 4b. "/7-shaped" 
spectra are characteristic of quasi-periodic orbits with 
initial conditions very close to the stable fixed points 
[37]. The existence of these abrupt peaks can be ex- 
plained theoretically. In fact such peaks are infinities 
[13,28]. These infinities have to be distinguished from 
local smooth maxima, which may be present in some 
spectra. The outermost black solid line shown in Fig. 
4b is the ZVC. Figure 4c depicts the S(g) spectrum for 
the 1:2 resonant orbit with initial conditions as in Fig. 
3a, while in Figure 4d we see the two islands of invari- 
ant curves in the (x,p x ) phase plane. In Figure 4e we 
present the S(g) spectrum for the 2:3 resonant orbit 
of Fig. 3a. In this case, we observe three well defined 
U type structures, as much is the number of islands of 
invariant curves shown in Figure 4f. Moreover, looking 
more carefully we see that the left and the right spectra 
are exactly symmetrical about the g = axis, while the 
central spectrum lies on both sides of this axis. This 
indicates that two of the islands of invariant curves are 
symmetrical about the p x axis, while the third inter- 
sects the p x axis. In other words, we have the case of 
a quasi-periodic orbit, with a starting point on the p x 
axis. Figure 4g and 4h are similar to Fig. 4e and 4f 
but for the 3:4 resonant orbit of Fig. 3a producing four 
islands of invariant curves. 

Figure 5a illustrates the S(g) spectrum for the 3:5 
resonant orbit of Fig. 3a, while in Figure 5b the corre- 
sponding islands of invariant curves are presented. Sim- 
ilarly, in Figure 5c we see the S(g) spectrum for the 4:5 
resonant orbit of Fig. 3a and in Figure 5d the five is- 
lands of invariant curves in the (x,p x ) phase plane. Fur- 
thermore, Figure 5e depicts the S(g) dynamical spec- 
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Fig. 5 (a-h): (left pattern): The S(g) dynamical spectrum for the different types of resonant orbits and (right pattern) the 
corresponding sets of islands of invariant curves in the (x,p x ) phase plane. 
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Fig. 6 (a-h): (left pattern): The S(g) dynamical spectrum for the different chaotic components and the unified chaotic sea and 
(right pattern) the corresponding chaotic components and the unified chaotic sea in the (x,p x ) phase plane. 
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trum for the 5:7 resonant orbit of Fig. 3a, while in Fig- 
ure 5f we see the set of seven small islands of invariant 
curves. Finally, in Figure 5g one may observe the struc- 
ture of the S(g) spectrum for the 5:8 resonant orbit with 
initial conditions as in Fig. 3a. The corresponding set 
of the eight small islands is given in Figure 5h. For all 
the S(g) spectrums presented in Figs. 4 an 5 the initial 
conditions are as in Fig. 3a, while the integration time 
equal to 5000 time units. From the above results, we 
may conclude that the S(g) spectrum has the ability 
to identify every type of resonant orbit, since it pro- 
duces always as many U type structures as much as 
the total number of islands of invariant curves in the 
(x,p x ) phase plane. Moreover, it provides fast and very 
reliable results not only for simple resonant 2D orbits, 
but also for complicated resonant orbits of higher mul- 
tiplicity which correspond to sets of multiple islands of 
invariant curves. 

Our next step, is to deploy the S(g) spectrum in 
order to identify the several chaotic components of the 
2D dynamical system. Figure 6a shows the S(g) spec- 
trum for the first chaotic component (CC-1). We ob- 
serve a complicated, highly asymmetric structure, with 
a lot of large, small and abrupt peaks, indicating chaotic 
motion. The first chaotic component (CC-1) inside the 
ZVC is shown in Figure 6b. Similarly, Figure 6c depicts 
the S(g) spectrum for the second chaotic component 
(CC-2), while CC-2 is presented in Figure 6d. More- 
over, Figure 6e depicts the S(g) dynamical spectrum 
for the third chaotic component (CC-3), while CC-3 is 
presented in Figure 6f. In Figure 6g one may observe 
the S(g) spectrum for the unified chaotic sea which is 
shown in Figure 6h. Here, we must point out that Figs. 
6a, 6c, 6e and 6g reveal the sensitivity of the S(g) spec- 
trum. It is obvious that the structure shown in Fig. 6g is 
quite different from those shown in Figs. 6a, 6c and 6e. 
In other words, the S(g) spectrum is sensitive enough 
allowing us to know whether an orbit corresponds to a 
chaotic sea or a chaotic component. The initial condi- 
tions for the S(g) spectra corresponding to the chaotic 
components and the chaotic sea, shown in Fig. 6 are as 
in Fig. 3b, while the integration time is equal to 5000 
time units. 

Our final step, is to verify if the S(g) spectrum can 
provide reliable results regarding the evolution of 2D 
sticky orbits. Figure 7a shows the S(g) spectrum of a 2:3 
resonant 2D orbit producing three islands of invariant 
curves. The initial conditions are xq = 0,ya = 0,p x o = 
14.6, while the value of p y o is always found from the 
energy integral (8). The values of all other parameters 
are as in Fig. 1. Here, we observe three well defined U 
type spectra, indicating regular motion. Note, that the 
number of spectra is equal to the number of islands. 



Figure 7b shows the S(g) spectrum of an orbit starting 
near the above regular orbit. The initial conditions are 
x Q = 0,y = 0,p x0 = 10.5. Here, we observe three dif- 
ferent spectra with a lot of large and small peaks. This 
indicates that we have a sticky orbit and the sticky 
region is composed of three sticky islands. The sticky 
period is about T = 740 time units. Figure 7c shows the 
spectrum of the same orbit, when T = 790 time units. 
Here, the three spectra have joined in order to produce 
a unified structure. This indicates, that after the sticky 
period the test particle has gone to the chaotic sea. In 
Figure 7d we can see the S(g) spectrum of the orbit 
when T = 5000 time units. Here, we observe a spec- 
trum of a highly chaotic orbit with a lot of small, large, 
asymmetric and abrupt peaks. Note that, the sticky 
period obtained from the S(g) spectrum is very close 
to the sticky period obtained by the 5(c) spectrum in 
Paper II (see for more details in Fig. 3 in Paper II). 

Let us now follow the time evolution of another 2D 
sticky orbit using the S(g) spectrum. Figure 8a shows 
the S(g) spectrum of a 3:4 resonant 2D orbit producing 
four islands of invariant curves. The initial conditions 
are x = 0.1, j/o = 0,p x o = 8.5, while the value of p y0 is 
always found always from the energy integral (8). The 
values of all other parameters are as in Fig. 1. Here, 
we observe four well defined U type spectra, indicat- 
ing regular motion. Note, that the number of spectra 
is once more equal to the total number of islands. Fig- 
ure 8b shows the S(g) spectrum of an orbit starting 
near the above regular orbit. The initial conditions are 
xq = 0.1,yo = OtPxQ = 7.45. Here, we observe four dif- 
ferent spectra with a lot of large and small peaks. This 
indicates that we have a sticky orbit and the sticky 
region is composed of four sticky islands. The sticky 
period is about T = 950 time units. Figure 8c shows 
the spectrum of the same orbit, when T = 1000 time 
units. Here, the four spectra have joined in order to pro- 
duce a unified structure. This indicates, that after the 
sticky period the test particle has left completely from 
the sticky region and moved to the chaotic sea. In Fig- 
ure 8d we can see the S(g) spectrum of the same orbit 
when T = 5000 time units. Here, we observe a highly 
chaotic spectrum with a lot of small, large, asymmet- 
ric and abrupt peaks. Therefore, it is evident that the 
S(g) spectrum can be deployed in order to calculate the 
sticky period of a 2D sticky orbit and also to follow its 
time evolution towards the chaotic sea. 

So far, we have seen that the S(g) dynamical spec- 
trum can provide very reliable results when we study 
the character of the orbits in a 2D Hamiltonian sys- 
tem. It has the ability to identify complicated resonant 
orbits, different chaotic components and also the evolu- 
tion of sticky orbits. However, all the above results ob- 
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Fig. 7 (a-d): {a-upper left): The S(g) dynamical spectrum of a 2:3 resonant quasi-periodic 2D orbit and (b-d): the time evolution 
of the S(g) spectrum for a nearby sticky orbit. See the text for details. 
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Fig. 8 (a-d): (a,-upper left): The S(g) dynamical spectrum of a 3:4 resonant quasi-periodic 2D orbit and (b-d): the time evolution 
of the S(g) spectrum for a nearby sticky orbit. See the text for details. 
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Fig. 9 The structure of the (x,p x ) phase plane reproduced 
using the S(g) spectrum. Inside the green regions the motion is 
ordered, while inside the red regions the motion is chaotic. The 
outermost black solid line is the Zero Velocity Curve (ZVC). 
(Color figure online) 

tained from the S(g) are qualitative. Therefore, we must 
ispect each time the shape of the spectrum's structure 
by eye, in order to characterize a 2D orbit. Obviously, 
this is not very practical when someone wants to check 
a large number of orbits, so as to form an idea about 
the global structure of the phase plane of a dynamical 
system. Thus, we need to establish a criterion, in order 
to quantify the results given by the S(g) spectrum. This 
criterion can be derived by looking the structure of the 
spectra shown in Figs. 4, 5, 6, 7 and 8. One may ob- 
serve, that when the orbit is regular the maximum value 
of S(g), that is S(g) max , is high (S(g) max > 7). On the 
contrary, when the orbit is chaotic S(g) max is signifi- 
cantly smaller {S{g) max < 3). By testing a large num- 
ber of orbits (approximately 1000), we conclude that 
when S(g) max > 2.5 the orbit is regular, while when 
S{g)max < 2.5 the orbit is chaotic. By using the above 
criterion for initial conditions {xo,p x $) on the entire 
phase plane of Fig. 1 and giving to each initial condi- 
tion a color according to the value of the S(g) max , we 
can have a clear picture of the regions where chaotic or 
ordered motion occurs. The outcome of this procedure 
for the 2D dynamical system (8), using a dense grid 
of initial conditions on the phase plane, is presented 
in Figure 9. Thus, in Fig. 9 we clearly distinguish be- 
tween green regions, where the motion is ordered and 
red regions, where it is chaotic. The outermost black 
solid line is the Zero Velocity Curve (ZVC). It is worth 
mentioning that in Fig. 9 we can observe small islands 
of stability inside the large unified chaotic sea, which 
are not visible in the detailed phase plane of Fig. 1. Al- 
though each orbit in Fig. 9 was computed for only 1000 
time units, this time interval was sufficient for the clear 
revelation of small ordered regions inside the chaotic 



sea. For a grid of about 6 x 10 4 equally spaced ini- 
tial conditions (%o->Pxo)i we need about 7 hours of CPU 
time on a Pentium Dual-Core 2.2GHz PC, in order to 
construct the grid presented in Fig. 9. 

3 Application of the S(w) dynamical spectrum 
in the 3D Hamiltonian system 

In this Section, we shall present results regarding the 
character of motion in the 3D dynamical system. For 
this purpose we use the 3D Hamiltonian (6) and take 
h = h 2 , that is the value of energy is equal to that of the 
2D system. Now, we can use the results obtained from 
the 2D system, in order to investigate the motion in 
the 3D system. What we do is to consider orbits with 
a starting a point on the (x,p x ) phase plane with an 
additional value of z = zq and follow the evolution of 
the 3D orbits. 

Figure lOa-h presents eight typical resonant periodic 
orbits in the 3D system. In all resonant orbits the values 
of the initial conditions (xo,p x o) are the same as in the 
corresponding 2D orbits, discussed in Fig. 2a-h. For all 
3D orbits y = p z0 = 0, the value of p y o is always found 
from the energy integral (6), while the initial value of zrj 
is equal to 0.005. Note, that all orbits shown in Fig. 10a- 
h stay very close to the galactic plane, due to the small 
initial value of zq. Nevertheless, we observe that even for 
a small value of Zq the large variety of resonant orbits we 
studied in the 2D system in the previous Section, still 
exists also in the 3D dynamical system. The integration 
time for all 3D orbits shown in Fig. lOa-h is 100 time 
units. 

A basic problem in Hamiltonian systems of three 
degrees of freedom is the visualization of the 4D sur- 
faces of section. Let us assume the phase space of an 
autonomous Hamiltonian system, that has 6 dimen- 
sions, e.g. in Cartesian coordinates, (x,y, z,x,y, z). For 
a given value of the Hamiltonian a trajectory lies on 
a 5D manifold. In this manifold the surface of section 
is 4D. This does not allow us to visualize it directly. 
Several methods have been applied to overcome this 
problem in the past and we summarize them below. 

The structure of the 4D space phase space was ex- 
amined for the first time in the pioneer work of Froeschle 
[14,15]. In that work, he used stereoscopic views and the 
method of slices in order to understand and interpret 
the structure of the tori, that appeared at the neigh- 
borhood of stable periodic orbits. Similar methods have 
also been applied for studying the 3D projections of in- 
variant tori in the 4D surface of section or in the phase 
space of a 4D symplectic map [8,24,25]. The 2D pro- 
jections of such invariant tori have been examined on 
various 2D planes in detail [33,34]. In the present work, 




Fig. 10 (a-h): Typical resonant periodic orbits in the 3D system. The initial conditions for all the orbits and additional details 
are provided in the text. 
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we shall use the method proposed by Patsis & Zachi- 
las [27]. In this method, we first define the surface of 
section that we will use, e.g. y = with p y > 0. Then, 
we select a 3D subspace projection e.g. (x,p x , z) of the 
(x,p x , z,p z ) 4D phase space. If the 3D projection is a 
well formed 3D invariant torus then the initial condi- 
tions correspond to a regular orbit. On the other hand, 
if instead of a well defined torus, we observe a vague 
cloudy structure in the (x,p x ,z) projection, then the 
corresponding 3D orbit is chaotic. 

Another interesting approach is the method intro- 
duced by Pfenniger [29] This method has been used 
by several researchers in the past [4,32]. We take sec- 
tions in the plane y = 0, p y > of 3D orbits, whose 
initial conditions differ from the plane parent periodic 
orbits only by the z component. The set of the result- 
ing 4-dimensional points in (x,p x , z,p z ) phase space are 
projected on the (z,p z ) plane. If the projected points 
lie on a well defined curve, lets call it an "invariant 
curve", then the motion is regular, while if not, the 
motion is chaotic. The projected points on the (z,p z ) 
plane show nearly invariant curves around the periodic 
points at z — 0, p z — 0, as long as the coupling is weak. 
When the coupling is stronger the corresponding pro- 
jections in (z,p z ) plane display increasing departure of 
the plane periodic point, up to making a direct orbit 
a retrograde one and vise versa. Here, we must define 
what one means by direct and retrograde 3D orbit. If 
consequents in the (z,p z ) section of the 3D orbit drop 
in one of the two domains of the corresponding section 
of 2D orbits at the same value of the Jacobi integral 
Ej, we can distinguish between direct and retrograde 
motion. Orbits which visit both domains are intermit- 
tently direct or retrograde. 

For the study of 3D orbits, we introduced and used 
by Zotos [40] (hereafter Paper III) the S(w) spectrum. 
Here we must remind to the reader that the S(w) spec- 
trum is the distribution function of the w parameter 
AN(w) 



S(w) 



(12) 



NAw 

where AN(w) are the number of the parameters w in 
the interval (w,w + Aw), after N iterations. The pa- 
rameter w, is defined as 

(Xi -Pxi)~ {Zi -Pzi) 



Wi = 



Pyi 



(13) 



where (xj, Zi,p X i,p V i,p Z i) are the successive values of 
the (x, z,p X7 p y ,p z ) elements of the 3D orbit. By defi- 
nition, the w parameter is based on a combination of 
coordinates and velocities of the 3D orbit. It would be 
very interesting to test the effectiveness of the S(w) 
spectrum in the present dynamical system and find out 
if it has the ability to identify all the 3D resonant or- 
bits. Moreover, it would be extremely beneficial, if it 



could distinguish the different chaotic components and 
also if it is sensitive enough in order to provide reliable 
results regarding the evolution of the 3D sticky orbits. 

We shall first apply the S(w) dynamical spectrum, 
in order to identify the different types of 3D resonant 
orbits. The left side of Figure lla-h shows the S(w) 
spectrums for the 3D resonant orbits, while the right 
side of the same figure presents the corresponding 3D 
invariant tori in the (x,p X7 z) phase subspace. Figure 
11a shows the S(w) spectrum for the 1:1 resonant 3D 
orbit with initial conditions (x ,p x o) as in Fig. 3a. As 
expected, one can observe a well defined U type struc- 
ture, indicating regular motion. The corresponding 3D 
invariant torus is shown in Figure lib. Fig. lib depicts 
a well formed double-ring shaped 3D invariant torus. 
Figure 11c depicts the S(w) spectrum for the 1:2 reso- 
nant 3D orbit with initial conditions (xo,p x o) as in Fig. 
3a, while in Figure lid we see the two 3D islands of in- 
variant tori in the (x,p x ,z) phase subspace. In Figure 
lie we present the S(w) spectrum for the 2:3 resonant 
3D orbit with initial conditions (xq,p x q) as in Fig. 3a. 
In this case, we observe three well defined U type struc- 
tures, as much is the number of 3D islands of invariant 
tori shown in Figure llf. Figure llg and llh are similar 
to Fig. lie and llf but for the 3:4 resonant 3D orbit 
with initial conditions (xo,p x o) as in Fig. 3a produc- 
ing four 3D islands of invariant tori. Similarly, Figure 
12a illustrates the S(w) spectrum for the 3:5 resonant 
3D orbit with initial conditions (xo,p x o) as in Fig. 3a, 
while in Figure 12b the corresponding five 3D islands 
of invariant tori are presented. In Figure 12c we see 
the S(w) spectrum for the 4:5 resonant 3D orbit with 
initial conditions (x 0} p xn ) as in Fig. 3a and in Figure 
12d the five 3D islands of invariant tori in the (x,p x , z) 
phase subspace. Furthermore, Figure 12e depicts the 
S(w) dynamical spectrum for the 5:7 resonant 3D orbit 
with initial conditions (xo,p x o) as in Fig. 3a, while in 
Figure 12f we see the set of seven small 3D islands of 
invariant tori. Finally, in Figure 12g one may observe 
the structure of the S(w) spectrum for the 5:8 resonant 
3D orbit with initial conditions (x ,p x o) as in Fig. 3a. 
The corresponding set of eight 3D invariant tori is given 
in Figure 12h. 

For all the S(w) spectrums presented in Figs. 11 
an 12 the initial conditions (x ,p x o) are as in Fig. 3a, 
Uo = Pzo = 0, the value of p y0 was found always from 
the energy integral (6), while the initial value of z is 
0.005. The integration time for the S(w) spectra shown 
in Figs. 11 and 12 is 5000 time units, while for the 3D 
invariant tori is equal to 2 x 10 4 time units. From the 
above results, we may conclude that the S(w) spectrum 
has the ability to identify every type of resonant 3D or- 
bit, since it produces always as many U type structures 
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Fig. 11 (a-h): (left pattern): The S(w) dynamical spectrum for different types of resonant 3D orbits and (right pattern) the 
corresponding 3D islands of invariant tori in the (x,p x ,z) phase subspace. 
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Fig. 13 (a-h): (left pattern): The S(w) dynamical spectrum for the different chaotic components and the chaotic sea and (right 
pattern) the corresponding (x,p x ,z) projections on the phase subspace. 
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as much as the total number of 3D islands of invariant 
tori in the (x,p x ,z) phase subspace. Moreover, it pro- 
vides fast and very reliable results not only for simple 
resonant 3D orbits, but also for complicated resonant 
orbits of higher multiplicity which correspond to sets of 
multiple 3D islands of invariant tori. 

In Paper I, we have shown that the different chaotic 
components in the 3D dynamical system do not merge 
in order to form a unified chaotic component but in- 
stead they continue to exist separately for time intervals 
much larger than the age of the Universe. Therefore, 
our next step, would be to deploy the S(w) spectrum 
in order to identify the several chaotic components of 
the 3D dynamical system. Figure 13a shows the S(w) 
spectrum for the first chaotic component (CC-1). We 
observe a complicated, highly asymmetric structure, 
with a lot of large, small and abrupt peaks, indicating 
chaotic motion. The 3D toroidal structure of the first 
chaotic component (CC-1) is shown in Figure 13b. Sim- 
ilarly, Figure 13c depicts the S(w) spectrum for the sec- 
ond chaotic component (CC-2), while the 3D toroidal 
structure of CC-2 is presented in Figure 13d. More- 
over, Figure 13e depicts the S(w) dynamical spectrum 
for the third chaotic component (CC-3), while the 3D 
toroidal structure of CC-3 is presented in Figure 13f. 
In Figure 13g one may observe the S(w) spectrum for 
the unified chaotic sea. Things are quite different in 
the (x,p x ,z) projection shown in Figure 13h. Here, in- 
stead of a well defined toroidal structure, we observe 
a vague cloudy structure of random 3D points. This 
cloudy structure indicates that the corresponding 3D 
orbit is highly chaotic. Here, we must point out that 
Figs. 13a, 13c, 13e and 13g reveal the sensitivity of the 
S(w) spectrum. It is obvious that the structure of the 
S(w) spectrum shown in Fig. 13g is quite different from 
those shown in Figs. 13a, 13c and 13e. In other words, 
the S(w) spectrum is sensitive enough allowing us to 
know whether an 3D orbit corresponds to a chaotic sea 
or a chaotic component. The initial conditions (xo,p x o) 
for the S(w) spectra corresponding to the chaotic com- 
ponents and the chaotic sea, shown in Fig. 13 are as in 
Fig. 3b, while the initial value of Zq is 0.01. The inte- 
gration time for the S(w) spectra shown in Fig. 13 is 
5000 time units, while for the 3D (x,p x ,z) projections 
is equal to 5 x 10 4 time units. 

Our final step, is to verify if the S(w) spectrum can 
provide reliable results regarding the evolution of 3D 
sticky orbits. Extensive numerical calculations indicate 
that unfortunately, the S(w) dynamical spectrum is 
not sensitive enough in order to provide reliable results 
about the evolution of 3D sticky orbits. This fact led us 
in Paper I, to introduce the S(k) spectrum for the study 
of the 3D sticky motion. The S(k) dynamical spectrum 
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Fig. 14 A grid of initial conditions (xo,p x o) when Zo = 0.007 
constructed using the S(w) spectrum. Inside the green regions 
the 3D motion is ordered, while inside the red regions the 3D 
motion is chaotic. The outermost black solid line defines the 
equation / (x,p x ; Zo) = h. (Color figure online) 

is a modified version of the S(w) spectrum and it was 
designed especially for the exploration of the 3D sticky 
motion. Chaotic orbits can be very important in the 
self-consistency of a dynamical model, especially those 
being sticky to stable resonant structures for long time 
intervals (i.e t > 100 time units). These sticky orbits 
appear at the boundaries of transition from an island 
of stability to the surrounding chaotic zone. 

So far, we have seen that the S(w) dynamical spec- 
trum can provide very reliable results when we study 
the character of the orbits in a 3D Hamiltonian sys- 
tem. It has the ability to identify complicated resonant 
3D orbits and different chaotic components but it can 
not be applied in order to study the evolution of 3D 
sticky orbits. However, all the above results obtained 
from the S(w) are qualitative. Therefore, we must check 
each time the shape of the spectrum's structure by eye, 
in order to characterize a 3D orbit. Thus, as we did 
in the previous section for the S(g) spectrum, we need 
to establish a numerical criterion, in order to quantify 
the results provided by the S(w) spectrum. Once more, 
this criterion can be derived by looking the structure 
of the spectra shown in Figs. 11, 12 and 13. One may 
observe, that when the orbit is regular the maximum 
value of S(w), that is S(w) max , is high (S(w) max > 7). 
On the contrary, when the orbit is chaotic S(w) max is 
significantly smaller (S(w) ma x < 3). By testing a large 
number of 3D orbits (approximately 1000), we conclude 
that when S(w) max > 3.5 the corresponding 3D orbit 
is regular, while when S{w) max < 3.5 the 3D orbit is 
chaotic. By using the above criterion for 3D orbits with 
initial conditions (xo-,p x o) on the entire phase plane of 
Fig. 1 with an initial value z$ = 0.007 and giving to 
each initial condition a color according to the value of 
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the S(w) ma x, we can have a clear picture of the regions 
where chaotic or ordered 3D motion occurs. The result 
of this procedure for the 3D dynamical system (6), using 
a dense grid of initial conditions (x ,p xa ,z = 0.007), 
is presented in Figure 14. Thus, in Fig. 14 we clearly 
distinguish between green regions, where the motion is 
ordered and red regions, where it is chaotic. The outer- 
most black solid line defines the equation 

f{x,p x ;zo) = ^p 2 x + V tot (x; z ) = h. (14) 

Although each 3D orbit in Fig. 14 was computed for 
only 1000 time units, this time interval was sufficient 
enough for the clear revelation of the dynamical struc- 
ture of the 3D phase subspace, when z = 0.007. The 
structure of the grid shown in Fig. 14 has a major dif- 
ference from the grid shown in Fig. 9 corresponding 
to the 2D model. We observe in Fig. 14 that the two 
large central regions of stability shown in Fig. 9, which 
correspond to the 1:1 resonant orbits now have disap- 
peared in the 3D system. This was expected because 
as we have shown in Paper I, 3D orbits with initial 
conditions (xq,p x q) inside these islands remain regular 
only when zq < 0.0065. For all permissible larger ini- 
tial values of zo the 3D orbits become chaotic. In Fig. 
14 the initial value of zo is 0.007. Thus, we verified not 
only the effectiveness but also the reliability of the S(w) 
spectrum. For a grid of about 6 x 10 4 equally spaced 
initial conditions (xo,p x o, zq = 0.007), we need about 
7.5 hours of CPU time on a Pentium Dual-Core 2.2GHz 
PC, in order to construct the grid presented in Fig. 14. 

4 Discussion and conclusions 

In the present paper, we applied two different types of 
dynamical spectra in order to study the properties of 
motion in Hamiltonian systems of two (2D) and three 
(3D) degrees of freedom. Our main objective, was to 
test, verify and prove the effectiveness and also the re- 
liability of these new indicators. For our tests, we chose 
the dynamical system described by the Hamiltonian 
(6). Our choice can be justified due to the fact that 
this potential displays a large variety of resonant orbits, 
different chaotic components and also several sticky re- 
gions in both 2D and 3D systems. Therefore, it provide 
us an excellent opportunity to test and also reveal the 
great efficiency of our new dynamical parameters. 

For the Hamiltonian system of two degrees of free- 
dom (2D), we used the S(g) dynamical spectrum. The 
main conclusions of our numerical calculations in the 
2D system can be summarized as following: 

1. The S(g) dynamical spectrum has the ability to 
identify any type of complicated resonant 2D orbits, as 



it produces as many [/-shaped structures as the total 
number of the islands of invariant curves in the (x,p x ) 
phase plane. Moreover, the numerical outcomes indi- 
cate that the S(g) spectrum can identify islandic mo- 
tion corresponding to higher multiplicity orbits which 
produce sets of several islands of invariant curves in the 
phase plane, no matter how tiny these islands are. 

2. In a 2D dynamical system where different chaotic 
components do exist, the S(g) dynamical spectrum is 
an excellent indicator in order to distinguish between 
chaotic components and unified chaotic domains. In 
particular, the structure of the S(g) spectrum corre- 
sponding to a chaotic component has many differences 
from the structure of the S(g) spectrum of a highly 
chaotic orbit. Therefore, the S(g) spectrum is sensitive 
enough regarding the chaoticity of a 2D orbit. 

3. One more advantage of the S(g) spectrum, is that 
it can be deployed in order to calculate the sticky period 
of a 2D sticky orbit and also to follow its time evolution 
toward the chaotic sea. In the transition regions from 
ordered to chaotic motion, the pattern of the S(g) spec- 
trum has a small number of abrupt peaks near stable 
periodic orbits. However, it develops a large number 
of peaks near the KAM boundaries separating them 
from the surrounding unified chaotic domain. The S(g) 
spectra of chaotic orbits near, but outside the KAM 
boundaries have also a large number of peaks similar 
to the spectra of orbits near, but inside the same KAM 
boundaries (see Figs. 7 and 8). 

In the case of the Hamiltonian system of three de- 
grees of freedom (3D), we applied the S(w) dynami- 
cal spectrum. Our main conclusions of the numerical 
experiments in the 3D system can be summarized as 
following: 

1. The S(w) dynamical spectrum has the ability 
to identify any type of complicated resonant 3D or- 
bits, as it produces as many [/-shaped structures as 
the total number of the islands of 3D invariant tori in 
the (x,p x , z) projection of the phase subspace. Further- 
more, our numerical outcomes indicate that the S(w) 
spectrum can identify 3D islandic motion correspond- 
ing to higher multiplicity orbits which produce sets of 
islands of 3D invariant tori in the phase subspace, no 
matter how tiny these 3D islands are. 

2. In a 3D dynamical system where different chaotic 
components still do exist, the S(w) dynamical spectrum 
is an appropriate indicator in order to distinguish be- 
tween chaotic components and unified chaotic domains. 
In particular, the structure of the S(w) spectrum corre- 
sponding to a chaotic component has many differences 
from the structure of the S(w) spectrum of a highly 
chaotic 3D orbit. Therefore, the S(w) spectrum is sen- 
sitive enough regarding the chaoticity of a 3D orbit. 
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Fig. 15 (a-d): (a-upper left): The S(ct) dynamical spectrum, (b-upper right): the S(r) dynamical spectrum, (c-lower left): The 
S(c) dynamical spectrum and (d-lower right): The S(g) dynamical spectrum for a 2D orbit. Details are provided in the text. 



3. Our numerical calculations indicate that unfor- 
tunately, the S(w) dynamical spectrum is not sensitive 
enough in order to provide reliable results about the 
evolution of 3D sticky orbits. This fact led us in Pa- 
per I, to introduce the S(k) spectrum for the study of 
the 3D sticky motion. Here we must point out, that the 
S(k) spectrum is a modified version of the S(w) spec- 
trum and it was constructed especially for the study of 
3D sticky motion. 

Thus, taking into account all the above findings, we 
conclude that the S(g) and S(w) spectrums proved to 
be very reliable and fast indicators which can help us to 
identify the character of orbits in Hamiltonian systems. 
Both spectra need about 5000 time units of integration 
time, in order to reveal the true nature of an orbit. 
Moreover, the definitions of these spectra are based on 
a complicated combination of the coordinates and the 
momenta (orbital elements), yielding more information 
about the structure of the orbit. Both spectrum defini- 
tions have the value of p y in their denominators. As the 
value of p y is always positive (p y > 0), while x, z, p x 
and p z can obtain positive and negative values, we can 
check the symmetry and locate the position of the cor- 
responding 2D islands of invariant curves in the (x,p x ) 
phase plane or the 3D invariant tori in the (x,p x ,z) 



phase subspace. In addition, the fact that p y is always 
positive is a great advantage because we avoid cases 
such as p y = which would cause problems in the def- 
initions of the g and w parameters. 

At this point, it would be of particular interest to 
compare our new spectrums with previously used spec- 
tra. The S{a) spectrum [37] has very limited ability to 
monitor the evolution of sticky motion. Furthermore, 
the S(a) spectrum needs the simultaneously computa- 
tion of two neighboring orbits, while our new spectra 
do not have this requirement. Another weakness of the 
S(a) spectrum is that it cannot distinguish complicated 
resonant orbits producing multiple sets of islands of in- 
variant curves. This is very important because (i) the 
detection of sets of small islands of invariant curves is 
significant in order to search for higher order resonances 
and (ii) there are cases where we have sets of small is- 
lands of invariant curves embedded in a large and uni- 
fied chaotic area. In these cases we must be able to 
detect the small islands of stability which exist as land- 
marks of ordered motion. Moreover, the S (c) spectrum 
[3] in several cases produces less U shaped structures 
than the total number of the islands of invariant curves 
or the invariant tori. Note, the coupling of the third 
component z, carrying all the information about the 
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3D motion, is hidden in the definition of the 5(c) spec- 
trum. Thus, the 5(c) spectrum provides implicit results 
regarding the 3D orbits. On the other hand, using the 
new S(w) spectrum the outcomes are explicit, since the 
coupling of the third component z is located directly in 
its definition. Therefore, by introducing the new S(g) 
and S(w) dynamical spectra, we have managed to over- 
come all the above drawbacks of the previously used 
spectra. 

It would be very illuminating for the reader, if we 
display an example comparing the structure of differ- 
ent kinds of dynamical spectra. We will present the 
dynamical spectra using chronological order. In Figure 
15 a-d we observe the structure of four different spec- 
tra corresponding to a 2D orbit with initial conditions 
x a = 0, yo = and p x o = 12.5, while the initial value 
of p y is found from the energy integral (8). In Figure 
15a we see the 5(a) spectrum [37]. It is clear that the 
5(a) spectrum is quite simple and almost symmetrical 
with respect to the a = axis. It tends to become U- 
shaped, but with a lot of small peaks. Moreover, since 
this orbit produces three islands of invariant curves in 
the (x,p x ) phase plane one should also expect three sep- 
arate [/-shaped structures. On the contrary in Fig. 15a 
we observe that the three spectra are joined together. 
Figure 15b displays the output of the 5(r) dynamical 
spectrum [20]. Here, we see that the spectrum again 
is almost symmetrical with respect to the r = axis. 
However, with a more closer look one may observe that 
the [/-shaped structures shown in Fig. 15b are not very 
smooth. Furthermore, in some cases the total number 
of [/-type structures is smaller than the correspond- 
ing number of islands of invariant curves, due to the 
fact that symmetric islands produce identical spectra. 
In order to overcome this drawback, a new type of dy- 
namical spectrum was introduced in [3]. In Figure 15c 
we present the structure of the 5(c) spectrum. In this 
case, the [/-type structures are very smooth but they 
are completely asymmetric with respect to the c = 
axis. Indeed, the main disadvantage of the 5(c) spec- 
trum is its inadequateness to reveal the relative posi- 
tions of the corresponding islands of invariant curves. 
In addition, it does not provide reliable results when 
the orbit is complicated and corresponds to a set of 
several small islands of invariant curves. The structure 
of the new S(g) spectrum is shown in Figure 15d. It is 
obvious, that here everything are as they should have 
been. The [/-type structures are separate, very smooth 
and completely symmetrical with respect to the g = 
axis. Thus, one may conclude that the S(g) spectrum 
is a powerful tool with impressive performance, when it 
is applied to detect resonant orbits of higher multiplic- 



ity. The integration type for all the spectra presented 
in Fig. 15a-d is 5000 time units. 

There are also other dynamical indicators, such as 
the LCE, FLI [17], SALI [35] and GALI [36], which can 
be used in order to characterize an orbit. The LCE is 
a well known dynamical indicator, but it requires vast 
time intervals of numerical integration of order of 10 5 
time units in order to provide reliable and conclusive re- 
sults about the nature of an orbit. On the contrary, FLI, 
SALI and GALI need at most 500 to 1000 time units of 
numerical integration time, in order to reveal the true 
nature of an orbit in Hamiltonian systems. The above 
methods can not be used in order to identify compli- 
cated resonant orbits of higher multiplicity. Moreover, 
FLI and GALI are not very sensitive in the case of sticky 
motion. However the FLI method is more sensitive in 
the case of thin chaotic layers. Furthermore, the SALI 
method can provide reliable information about sticky 
orbits, using shorter integration time than the S(g) or 
S(w) spectrums. The main disadvantage of FLI, SALI 
and GALI is that these indicators need the simultane- 
ously integration of several sets of variational equations 
and deviations vectors, while the S(g) and S(w) dy- 
namical spectra need only the integration of the basic 
equations of motion. 

All the different spectral definitions used in previous 
studies could obtain only qualitative results regarding 
the nature of an orbit. Thus, every time we needed to in- 
spect the shape of the spectrum's structure by eye in or- 
der to characterize an orbit. In the present research we 
have established numerical criteria in order to quantify 
the results obtained from the S(g) and S(w) dynamical 
spectra. It was shown in Figs. 9 and 14 that these nu- 
merical criteria are valid and therefore, using the S(g) 
and S(w) dynamical spectra we can reveal the dynam- 
ical structure of the phase plane through the construc- 
tion of a dense grid. In addition, the S(k) spectrum is 
a very effective tool for the study of 3D sticky motion. 
Thus, it is the author's humble opinion, that by using 
the triplet of (S(g), S(w), S(k)) dynamical spectra, one 
can study and obtain very fast and reliable results, re- 
garding the structure of Hamiltonian systems of two 
or three degrees of freedom. Since these new dynami- 
cal spectra have been proved extremely accurate and 
reliable, they may not be amenable to further improve- 
ments and thus, we may have marked the end of an 
era for the dynamical spectra of orbits in Hamiltonian 
systems. 
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